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Sum rule tests are performed on the spectral data for (i) flavor ud vector-current- 
ly induced hadronic r decays and (ii) e + e~ hadroproduction, in the region below 
s ~ 3 — 4 GeV 2 , where discrepancies exist between the isospin-breaking-corrected 
charged and neutral current 1 = 1 spectral functions. The r data is found to be 
compatible with expectations based on high-scale a s (Mz) determinations, while the 
Q | electroproduction data displays two problems. The results favor determinations of 
\Q ' the leading order hadronic contribution to (g — 2)^ which incorporate hadronic r 
decay data over those employing electroproduction data only, and hence a reduced 
■ discrepancy between experiment and the Standard Model prediction for (g — 2) M . 

> 

^ ' PACS numbers: 13.35.Bv,13.40.Gp,13.66.Bc,13.35.Dx 
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O ■ I- INTRODUCTION 

Ph! In the Standard Model (SM), the largest of the non-purely-leptonic contributions to 

the anomalous magnetic moment of the muon, a M = (g — 2)^/2, is that due to the leading 



order (LO) hadronic vacuum polarization (VP), [c^]^. a M is currently known to 0.5 
ppm [H, with a proposal to reduce this to 0.2 ppm in the near future 0- The 0.5 



ppm uncertainty represents < 1% of making an accurate determination of 

| crucial to the study of possible non-SM contributions to a^. 

\. a t^had * s re l & ted to a [e + e~ — > hadrons] by the dispersion representation 

lo «i A /(o) r J K (s 



= / ds^R(s), (1) 

where the form of K(s) is well-known, and R{s) is the "bare" e + e~ — > hadrons to 
e + e _ — > cross-section ratio. With recent electroproduction data, the uncertainty 
on [afj] had from Eq. Q is comparable to the experimental error on a M , and dominates 
the uncertainty in the SM prediction for 0, IE El HI- Since CVC relates the 1 = 1 
electromagnetic (EM) spectral function to the charged current vector spectral function 
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measured in r — > u T + non-strange hadrons, hadronic r decay data 0, 0] can, in 
principle, be used to improve the determination of 0, EH- At the < 1% level 

necessitated by the current experimental error, isospin-breaking (IB) corrections must be 
taken into account. 

IB corrections for the nn final state, whose contributions dominate [ajf°,, were stud- 
, |15J] and, together with kinematic IB corrections for the Air contribution, 
incorporated into the latest r-based [a^ad analyses 0,0]. Even after IB corrections, how- 
ever, the high-precision CMD-2 -kit EM data 13 lies ~ 5 - 10% below the IB-corrected 
t data for between ~ 0.85 and ~ 1 GeV 0, [lTf. The corresponding determina- 
tions of {a^ad differ by ~ 2a, the r-based result lying higher andproducing a SM 
prediction in significantly better agreement with experiment (USB HI- Recent KLOE 
e + e _ — > 7r + 7r~ cross-sections [l9[ yield an [a M ]^ compatible with CMD-2 Q, though the 
point-by-point agreement between the two data sets is not entirely satisfactory [7[ . 

In view of the unsettled experimental situation, we study sum rule constraints on 
weighted integrals of the 1 = 1 vector r decay and EM spectral functions. The weights, 
w(s), and upper integration limits, so, are chosen such that (i) each spectral integral has 
a reliable and well-converged OPE representation, and (ii) all relevant OPE input can be 
obtained from sources independent of the low scale EM and r data we seek to test. OPE 
uncertainties are minimized by working with s and w(s) for which the relevant OPE 
representation is dominated, essentially entirely, by its D = component, and hence 
determined, essentially entirely, by the single input parameter, a s (Mz), which can be 
taken from independent high-scale studies. The r decay based spectral integrals will be 
shown to be well reproduced by the corresponding OPE representations. Those based 
on EM data, in contrast, will be shown to lie consistently below the corresponding OPE 
values, for positive w(s), and to differ from them in their so dependence. Both features 
are as expected if the EM spectral data is too low in the disputed region. 



II. THE SUM RULE CONSTRAINTS 

We study sum rule constraints on the EM spectral function, Pem(s), and the sum of 
the spin J = and 1 components of the charged 1 = 1 vector current spectral function, 
Pvtd( s ) = Pv~ud( s ) + Pvm( s )- The former is related to R(s) by pem{s) = R{s)/12tt 2 , 
and to the bare e + e~ —> hadrons cross-sections, o~bare(s), by 

Pem(s) = s a bare (s) /167!- 3 a E Ai(0) 2 . (2) 

Defining R v . ud by R v . ud = T[t~ -> u T hadronsv ;ud (j)]/T[t~ -> u T e~u e {y)] and y T = 
s/m 2 ., PvXd( s ) is related to R v -ud by 

R v , ud = 12n 2 \V ud \ 2 S EW [ dy T (l-y T ) 2 \(l + 2y T )p { ^ d \s)-2y T p^ ud (s)] (3) 

Jo 1 J 



where V ud is the flavor ud CKM matrix element, and Sew is an electroweak correction [21]. 
Contributions to Py. ud (s) are of 0([m d — m u } 2 ), and hence numerically negligible, allowing 



3 



Pv-,ud ( s ) ^° be determined from the experimental decay distribution. 



A. Finite Energy Sum Rules 



For any correlator, U.(s), with no kinematic singularities, and any w(s) analytic in 
s\ < M with M > sq, analyticity implies the finite energy sum rule (FESR) relation 



where p(s) is the spectral function of n(s). In QCD, for very large so the OPE represen- 
tation can be employed on the RHS of Eq. (jlj). As so is decreased, this representation is 
expected to break down first near the timelike real s-axis [22] . A range of "intermediate" 
So will thus exist for which the OPE, though unreliable for general w(s), will remain 
valid for those w(s) satisfying w(s = s ) = 0. The corresponding FESR's are called 
"pinched" or pFESR's. For vector (V) and axial vector (A) correlators, and w(s) = s , 
OPE breakdown (duality violation) is significant at so ~ a few GeV 2 |23j|. Polynomials 
w(y) (with y = s/sq) having even a single zero at s = sq (y = 1), however, remove such 
violations for so greater than ~ 2 GeV 2 jHf , even for the flavor ud V-A correlator 0| • 

In interpreting FESR results, one should bear in mind that very strong correlations 
exist between spectral integrals corresponding to the same w(y), but different s . Such 
correlations are particularly strong when w(y) > over the relevant interval, < y < 1, 
and even more so when w(y) is monotonically decreasing. Similar strong correlations 
exist between spectral integrals corresponding to different w (y), but fixed so- Correlations 
among the corresponding OPE integrals are also very strong, especially when the OPE 
is dominated, as below, by a single (in this case, D = 0) contribution. 

We work, in what follows, with pinched polynomial weights, w(y) = 'Yl im c m y in - The 
pinching condition, w{l) = implies ^2i m c m = 0. For reasons explained below, w(y) 
is further restricted to be non- negative and monotonically decreasing on < y < 1. 
Since Sq is the only scale entering the RHS of Eq. (@J), it is obvious, on dimensional 
grounds, that integrated OPE contributions of dimension D = 2k + 2 scale as 1/sq, up 
to logarithms. For D > 2, such contributions are absent (up to corrections of 0(a s )) 
when C(£>_2)/2 = 0. The structure of the logarithmic integrals, ^j s | =so ds y k £n(Q 2 / fi 2 ) / Q D , 
responsible for the 0(a s ) corrections, is such that cancellations inherent in the pinching 
condition J2 m c m = lead to strong numerical suppressions of these corrections. D > 8 
contributions are typically assumed to be negligible, since the relevant condensate values 
are not known phenomenologically. The much stronger s dependence of such high D 
contributions allows this assumption to be tested explicitly. 

The reason for working with non- negative, monotonically decreasing w(y) is that the 
EM spectral data for the n + 7i~ and 7r + 7r~7r°7r° states, which dominate the EM-r discrep- 
ancy, lie uniformly below the IB-corrected r data. Non- negativity of w(y) then ensures 
that, for all sq, the normalization of the r-based spectral integrals will be too high if it is 
the EM data which is correct, while the normalization of the EM spectral integrals will 
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be too low if it is the r data which is correct. Since the y value for a given experimental 
bin decreases with increasing s > a monotonically decreasing w(y) similarly ensures that 
the slope with respect to Sq of the r spectral integrals will be too high if the EM data is 
correct, while the slope with respect to sq of the EM spectral integrals will be too low if 
the r data is correct. The slope constraint is particularly useful because the slope of the 
corresponding OPE integrals is very tightly constrained, and only very weakly dependent 
on the dominant OPE input parameter cx s (Mz)- 



B. OPE Input 

Re-writing the weighted pFESR OPE integrals of the relevant correlator, II, in terms 
of the Adler function D(Q 2 ) = — Q 2 dU(Q 2 )/dQ 2 , allows potentially large logs in the 
D = contribution to be summed point-by-point along the integration contour. This 
"contour-improved" (CIPT) prescription is known to improve the convergence of the 
integrated D = series |25| . The Adler function is given by 

[ D (Q 2 )] D=0 = ^E4 0) « fc , (5) 

k>0 

i2\ _ „, ir\2\ 



where a = a(Q ) = a s (Q )/n is the running coupling in the MS scheme, and C — 1, 

2/3 for the r, EM cases, respectively. For N f = 3, d^ = = 1, df = 1.63982 
(°) _ a Q7im lidl A „nr;^„+„ A®) A®) 



and fig = 6.37101 [26]. An estimate for d\ , d\ = 27 ± 16 12711 also exists, based on 



methods which (i) work well for the coefficients of the D = series j2Sj and (ii) produced, 
in advance of the actual calculation, an accurate prediction for the recently computed 
0(a 3 ) D = 2 coefficient of the (J) = (0 + 1) V+A correlator sum 0. 

The leading D = 2 contributions for the r correlator are 0(m^ d ) and numerically neg- 
ligible. For the EM correlator, up to tiny 0{m 2 Ld /m 2 ) corrections, the D = 2 contribu- 
tions are determined by a and the running MS strange mass rh s . At the scales employed 
here the integrated D = 2 contribution is small. The full expression for [U EM (Q 2 )] D=2 
may be found in Ref. 3(|. The D = 4 terms in the OPE of the EM and r correlators 



are determined by the RG invariant light quark, strange quark and gluon condensates, 
(££)rgi, (ss)rgi and (aG 2 )f/nT : up to numerically tiny O(mf) corrections. The expres- 
sions may be found in Refs. (31 , 32]. The integrated D = 4 contributions are again small 



at the scales employed. To reduce OPE uncertainties, we concentrate here on weights 
for which the integrated leading D = 6 contributions are absent. More extensive studies 
will be reported elsewhere [33]. We assume throughout that contributions with D > 8 
may be neglected, but check this assumption for self-consistency, as discussed above. 

For D = A input we use (i) (aG 2 ) = (0.009 ± 0.007) GeV 4 (from the recent re- 
analysis of charmonium sum rules [34]) and (ii) (2mg££) = —m 2 f 2 (the GMOR relation). 
(m s ss)nGi then follows from conventional ChPT quark mass ratios and the standard 
estimate (ss)rgi/ (££)rgi — 0.8±0.2. For the D = 0, 2 input, a and m s , we employ exact 
solutions based on the 4-loop-truncated f3 and 7 functions [35|, with initial conditions 

m s (2 GeV) = 95 ± 20 MeV (6) 
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a s {M z ) = 0.1200 ±0.0020. (7) 

Eq. (JOJ) reflects the range of results obtained in recent sum rule and N f = 2 
and 2 + 1 unquenched lattice studies [3^ . The Nf = 5 value in Eq. (JJJ) is run down to the 
Nf = 3 low-scale region using standard 4-loop running and matching 13611 . with the same 
matching scales as used in the recent EM sum rule studies of Refs. jy, |20( (HMNT). The 
input a s {Mz) in Eq. © differs from the PDG 2004 average for the following reasons. 
First, the PDG average includes hadronic r decay input, which must be excluded if we 
wish to test the r decay data. Second, the PDG average is strongl y af fected by the quoted 
low, small-error determination from heavy quarkonium decay [39]. The Quarkonium 
Working Group, however, has (i) strongly criticized the input to the low central value, (ii) 
argued that the quoted error is underestimated by a factor of 3 — 5 ,41], and (iii) concluded 
the method is not competitive with extractions based on perturbative treatments of high 
scale processes ji^ . Eq. (j7J) is obtained by removing lower scale determinations, including 
those based on heavy quarkonium and r decay, from the PDG average [43j]. 



C. Spectral Input 

The ALEPH and CLEO r decay distributions are in good agreement. For definiteness, 
we use the ALEPH results, for which the covariance matrix is publicly available. A 
small global rescaling accounts for minor changes in B e , B^, and the strange branching 
fraction, B s , since the original ALEPH publication PDG04 [4(j values are used for 
B e and B^, while the updated B s value incorporates (i) the new (2004) world averages for 
B[t~ ->■ v r K-^\ and B[t~ -> v T K~it^~7r~] Q, (ii) the new (2005) CLEO results for the 
branching fractions of four-particle modes with kaons (45[, and (iii) the higher precision 
K^2 value for the K pole contribution j4(j ■ The long-distance EM corrections determined 
in Ref. |l5| are also applied to the dominant r~ — ► 7r~7r° contribution ^| [i^ j. 

Detailed discussions and assessments of the EM hadroproduction data base, as of 
2002-2003, can be found in Refs. 0, H 0. The exhaustive compilation of Ref. 
provides useful information on the treatment of radiative and VP corrections for older 
experiments where such details are absent from the original publications. Detailed as- 
sessments of the needed residual VP corrections are also contained in Refs. 0, These 
corrections are computed using the most recent version of F. Jegerlehner's code |48[, gen- 
erously provided by its author. We also employ the following new, high-precision results, 
published subsequent to the analyses of Refs. 0, 0], and the compilation of Ref. |4^ |: 
(i) the final published version of the SND An cross-sections |49[ (with systematic er- 
rors significantly reduced over those of the earlier preprint version); (ii) the updated 
CMD-2 tt+tt-tt+tt- and tt°7, r]j cross-sections; and (iii) the BABAR 3vr 
and 7T + 7r~7r + 7r~ (53| cross-sections. Since the nn component of the EM-r discrepancy is 
driven by the CMD-2 data, with its very small 0.6% systematic error [16], we employ 
only CMD-2 data in the CMD-2 region. Where the existence of newer data permits, we 
exclude older data for which systematic errors are incompletely stated, or absent, in the 
original publications, and/or the residual radiative/ VP corrections to be applied are un- 
known. Fortunately, data with missing systematic errors for which no such replacement 
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is possible play only a small numerical role. We assign a guess of 20% for these errors 
in such cases. For the small number of remaining older experiments where the situation 
with regard to residual VP corrections is unclear, we apply no VP correction. In all such 
cases, however, (i) the corresponding contributions to the spectral integrals are small, 
and (ii) the neglected VP corrections are, in any case, much less than the quoted sys- 
tematic errors. The treatment of "missing mode" contributions, and the use of isospin 
relations for a number of small contributions where direct experimental determinations 
are absent, or have large errors, follow the treatments discussed in detail in Refs. |4j, |6(. 
More details on the treatment of the EM data will be provided elsewhere j33[. 



III. ANALYSIS AND RESULTS 

For reasons discussed above, we concentrate on pFESR's involving non-negative, 
monotonically decreasing w(y). The only such degree 1 weight is Wi(y) — 1 —y. Weights 
with a double zero at y = 1, which more strongly suppress OPE contributions from 
the vicinity of the timelike axis, should be even safer from the point of view of poten- 
tial duality violation. A convenient set of such "doubly-pinched" weights is the family, 
w N(y) = 1 — jtzt V + jfzi y N ■, N > 2. For a given wn, the only non-a s -suppressed ("un- 
suppressed" in what follows) D > 4 OPE contribution is that with D = 2N + 2. This 
contribution scales as 1/ Sq +1 relative to the leading integrated D = term. The strong 
s -dependence allows the neglect of D > 8 contributions (for wn>?) to be tested for self- 
consistency. We have also studied pFESR's based on a number of other weights. Since 
the results in all cases point to the same conclusion, and OPE uncertainties are reduced 
for weights having no unsuppressed D = 6 contribution, we focus on the pFESR's for two 
of the weights defined above, w± and wq. Other results will be presented elsewhere |33| . 

Combined OPE errors for the various pFESR's are obtained by adding in quadrature 
uncertainties associated with the OPE input parameters (D — 4 condensates, m s , and 
a s (Mz)) and the truncation/residual scale dependence of the integrated D = 0, 2 series. 
The latter are estimated to be twice the magnitude of the last term in the corresponding 
truncated series. The resulting OPE error estimate is somewhat more conservative than 
that employed by HMNT. 

Results for the EM case are presented in Figures [TJ |2l those for the r case in Figures E3 
0] The dashed lines represent the central OPE results, the solid lines the upper and lower 
edges of the OPE error bands. Because of the strong correlations, the OPE band is better 
thought of as a bundle of allowed parallel lines than as a generally allowed region. We 
see that, in the region 2 GeV 2 < so < both the magnitude and slope of the integrals 
over the r decay distribution are in good agreement with OPE expectations. In contrast, 
the EM spectral integrals are consistently low relative to OPE expectations (particularly 
for So greater than ~ 2.5 GeV 2 ) and have slopes with respect to So significantly lower 
than those of the OPE curves. Both features of the EM results are as expected if the r 
data is correct, and hence the EM data low, in the disputed region. 

The implications of the normalizations of the EM and r spectral integrals can be 
quantified by working out the a s {Mz) required to match the OPE and spectral sides of 
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TABLE I: The values of a s (Mz) obtained by fitting to the experimental EM and r spectral 
integrals for sq = rn%., with central values for the D = 2,4 OPE input 



Weight 


MM Z )} EM 


MM z )] r 


w 6 


n 1 1 qq+0.0030 
n 1 1 ^n+0.0022 

u.nou_ 00026 


n 191 o+0. 0027 
U.1Z1Z_q 0032 

fl 1 1 qc;+0.0020 

U.ll»d_Q 00 22 



a given pFESR. The reliability of the OPE, and hence of the extraction of a s (Mz), is 
optimized by choosing so as large as possible - in the case of hadronic r decay, so = m%. 
Results corresponding to central input for the small D = 2, 4 OPE contributions and 
the So = vn? T values of the spectral integrals are shown in Table |l] for both the EM and 
r cases. The agreement between the r-decay and independent high-scale determinations 
is excellent [54]. In contrast, the EM data corresponds to a s (Mz) ~ 2a lower than the 
high-scale determination. The y 6 term of we{y), in principle, produces an unsuppressed 
D = 14 OPE contribution scaling as 1/sq (1/sq relative to the leading D = term). 
Such a contribution, if present, would contaminate the extraction of a s (M z )- There 
is, however, no evidence for such a contribution, at a level which would impact our 
analysis, in the so dependence of either the EM or r u> 6 -weighted spectral integrals |5oT |. 
The excellent agreement between the a s [Mz) extracted using different doubly-pinched 
weights, with potential unsuppressed D > 6 contributions of different dimension, provides 
further evidence in support of the absence of such D > 6 contributions [5a, |57| . 

To quantify the disagreement between the EM OPE and experimental slope values, we 
work out the correlated errors for the slopes with respect to Sq of the OPE and spectral 
integrals. The correlations are such that the uncertainty on the OPE side is rather small. 
In particular, the slope is quite insensitive to a s [Mz)- These points are illustrated in 
Table ITT1 which shows the spectral integral and OPE slope values for the EM w\ and wq 
pFESR's. The OPE entries labelled "indep" are those obtained using the independent, 
high-scale fit value for a s (Mz)- Those labelled "fit", in contrast, correspond to the 
a s (M z ) values obtained by fitting to the relevant s = spectral integrals, as given in 
Table |U We see that, even if one were willing to tolerate the lower central a s {Mz) values 
implied by the EM spectral integrals, such a lowering of a s (Mz) would have negligible 
impact on the OPE vs. spectral integral slope discrepancy problem. 

TABLE II: Slopes wrt sq of the EM OPE and spectral integrals 



Weight 


Sexp 


a s {M z ) 


SoPE 


Wl 


.00872 ± .00026 


indep 


.00943 ± .00008 






fit 


.00934 ± .00008 


Wq 


.00762 ± .00017 


indep 


.00811 ± .00009 






fit 


.00805 ± .00009 
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IV. DISCUSSION AND CONCLUSIONS 

We have shown that weighted spectral integrals constructed using 1 = 1 hadronic 
t decay data are in good agreement with OPE expectations, while those involving EM 
data (i) require a value of a s {Mz) ~ 2<r below that given by high-scale determinations 
and (ii) correspond to a slope with respect to Sq in ~ 2.5cr disagreement with the OPE 
prediction. The slope problem, moreover, cannot be cured simply by adopting the lower 
a s (Mz) values, shown in Table HI which would bring the normalization of the OPE and 
spectral integrals into agreement for s — vn? T . The insensitivity of the slope to a s [Mz) 
also means that the agreement between the OPE expectation and the observed slope for 
the r decay spectral integrals represents a non-trivial test of the r data. 

One possibility is that the problems with the EM sum rules might be attributable 
to the presence of residual duality violation at the intermediate scales studied here; the 
success of the OPE in predicting both the slope and magnitude of the r-decay-based 
spectral integrals over the whole of the region 2 GeV 2 < sq < however, renders such 
an explanation highly implausible. The results thus point to the reliability of the r data, 
and to the likelihood of either (i) a problem with the experimental EM spectral distribu- 
tion, or (ii) the presence of as-yet-unidentified non-one-photon physics contributions in 
the experimental EM cross-sections. This in turn suggests that a M determinations which 
incorporate r decay data are to be favored over those employing EM spectral data only. 

While the disagreement between the EM and high-scale determinations of a s (Mz) is 
only ~ 2a, even with significantly lower high-scale input, e.g., the 2002 PDG average 
a s (Mz) = 0.1172 ± 0.0020 used by HMNT, the EM normalization, and even more so 
the EM slope, would still require, on average, upward fluctuations in pem(s). Since 
K(s)/s > 0, such fluctuations would typically also increase [a^]^ ad - Thus, even ignoring 
the slope problem and assessing the EM data as moderately consistent, within errors, 
with the OPE constraints, the fact that the spectral integrals lie consistently below 
the corresponding OPE constraint values, for any sensible input a s (Mz), points to the 
likelihood of a [o/Jw value higher than the current central EM-data-based value. 

Two further points are of relevance to assessing the implications of our results for the 
value of [o/JhocT Fi rs t> we find that, replacing the EM , 7r + 7r~7r°7r° and 

data with the equivalent r data resolves completely both the normalization and slope 
problems for the resulting modified "EM" spectral integrals j58[. Second, it is readily 
demonstrated that the pFESR's employed are sensitive to, not just the discrepancies 
in the An region, but also those in the 2tt region. (This is relevant since the [a^jf d 
integral is dominated by the 2n spectral contribution.) In fact, for the W\ pFESR, the 
shift in the EM spectral integral associated with the modification of the 2ix part of 
the EM spectral function represents 82% of the full shift at Sq = 2 GeV 2 and 32% at 
s = m 2 . The corresponding figures for the w 6 pFESR are 87% at s = 2 GeV 2 and 
45% at So = m^.. Thus, even though the pFESR's employed are relatively more sensitive 
to the 47r spectral contributions than is the [a^]f^ d integral, a clear sensitivity to the 
2tt component remains, making the constraints associated with these pFESR's highly 
relevant to the [o^], . problem. 
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In conclusion, all the sum rule tests performed favor the reliability of the r decay data, 
and point to problems with the EM data. We conclude that, at present, determinations 
of [dfj]^ d employing IB-corrected r decay data are more reliable than those based on EM 
data alone, and hence that there is no clear sign of a discrepancy between the current 
experimental value for a M and the SM prediction. 

NOTE ADDED: Subsequent to the submission of this paper, new results for the e + e~ — > 
7T + 7r _ cross-sections were released by the SND Collaboration L r 59]. As would be expected 
from the sum rule results above, the SND cross-sections are compatible with the r — > v t tttt 
data, but in significant disagreement with the KLOE tttt data, in the disputed region. 
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in excellent agreement with those obtained using the corresponding w$ sum rule. 
The ALEPH 9] and OPAL collaborations have, in fact, performed extractions of the 
effective D = 6, 8 condensate combinations for the ud V correlator by fitting to the sq = m% 
spectral integrals for a number of different "spectral weights" . The reader is cautioned that 
these extractions are performed assuming D > 8 contributions are absent, despite the fact 
that, for all but one of the weights employed, such contributions are, in principle, present 
in unsuppressed form. The analogous assumption for the ud V-A correlator has been tested 
explicitly and found to be incompatible with data 24 1, making it also suspect for the V 
case. One can nonetheless take the extracted values as representative of what might be 
expected for the size of such effects and investigate the impact on the sum rules used in 
this analysis. By construction, no unsuppressed D = 6 contributions are present in the 
sum rules considered. While the D = 8 term does not contribute to the w\ or w§ pFESR's, 
it does contribute to the u>3 pFESR. Using the larger of the ALEPH and OPAL D = 8 
fitted values, it is straightforward to demonstrate that, at sq = m^., the resulting D = 8 
contribution is at the level of 0.2% of the dominant D = term, producing a shift in the 
value of a s (Mz) extracted using the sum rule of < 0.0003, which, as claimed, is much 
smaller than the ~ .0020 error associated with data uncertainties. This observation simply 
confirms the conclusion obtained by studying the sq dependence of the various sum rules 
and the consistency of the a s {Mz) extractions from sum rules where non-a^-suppressed 
contributions with different D > 6 would be present |s3 |. 

The normalizations of the modified EM integrals are now in full agreement with the high- 
scale-input OPE constraint bands, for all so > 2 GeV 2 . The modified spectral integral 
slopes are also in good agreement with OPE expectations (the modified central values 
are 0.00947 for the w\ and 0.00807 for the w 6 pFESR, to be compared with the OPE 
predictions 0.00943 ± 0.0008 and 0.00811 ± 0.00009 already reported in Table HU 
M.N. Achasov et al. (The SND Collaboration), |arXiv:hep- ex/0506076 
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FIG. 1: EM OPE and spectral integrals for the weight w\ 
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